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Abstract 

One of the first and therefore most important theorems in perturbation theory 
claims that for an arbitrary self-adjoint operator A there exists a perturbation B 
of Hilbert-Schmidt class with arbitrary small operator norm, which destroys com¬ 
pletely the absolutely continuous (a.c.) spectrum of the initial operator A (von 
Neumann). However, if A is the discrete free 1-D Schrodinger operator and B is an 
arbitrary Jacobi matrix (of Hilbert-Schmidt class) the a.c. spectrum remains per¬ 
fectly the same, that is, the interval [—2,2]. Moreover, Killip and Simon described 
explicitly the spectral properties for such A + B. Jointly with Damanik they gen¬ 
eralized this result to the case of perturbations of periodic Jacobi matrices in the 
non-degenerated case. Recall that the spectrum of a periodic Jacobi matrix is a 
system of intervals of a very specific nature. Christiansen, Simon and Zinchenko 
posed in a review dedicated to F. Gesztesy (2013) the following question: “is there 
an extension of the Damanik-Killip-Simon theorem to the general finite system of 
intervals case?” In this paper we solve this problem completely. Our method deals 
with the Jacobi flow on GMP matrices. GMP matrices are probably a new ob¬ 
ject in the spectral theory. They form a certain Generalization of matrices related 
to the strong Moment Problem, the latter ones are a very close relative of Jacobi 
and CMV matrices. The Jacobi flow on them is also a probably new member of 
the rich family of integrable systems. Finally, related to Jacobi matrices of Killip- 
Simon class, analytic vector bundles and their curvature play a certain role in our 
construction and, at least on the level of ideology, this role is quite essential. 


1 Introduction 

1.1 Main result 

flj Von Neumann Theorem [42] states that for an arbitrary self-adjoint operator A, 
having a nontrivial absolutely continuous (a.c.) component of the spectrum, there 
exists a self-adjont perturbation dA of Hilbert-Schmidt class such that A + 5A has a 
pure point spectrum. Moreover, dA may have an arbitrary small operator norm. 

Therefore, the following result is already quite non-trivial. 

(2) Deift-Killip Theorem [11]. For a discrete one-dimensional Schrodinger operator with 
square summable potential, the absolutely continuous part of the spectrum is [—2,2]. 
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Thus, under special perturbations of Hilbert-Schmidt class (the square summable 
potential) the absolutely continuous spectrum of the free, discrete 1-D Schrodinger 
operator is perfectly preserved. It is totally surprising that one can find a complete 
explicit characterization of the spectral data if the perturbation is an arbitrary Jacobi 
matrix of Hilbert-Schmidt class. 

(3) Killip-Simon Theorem [20]. Let da be a probability measure on M with bounded but 
infinite support. As it is well known the orthonormal polynomials Pn{x) with respect 
to this measure obey a three-term recurrence relation 

xPn{x) = a{n)Pn-i{x) + b{n)Pn{x) -h a{n + l)Pn+i{x), a{n) > 0. 

The following are equivalent: 

(op) En>i Hn) - Ip < oo and < oo- 

(sp) The measure da is supported on [—2, 2] U X, and moreover 

log a'{x)\^/4 — x'^dx -b ^ \Jx\ — 4, < oo. (1.1) 

Xk&X 



Remark 1.1. Of course the (op)-condition means that the Jacobi matrix 




6(0) a(l) 
a(l) 6(1) a(2) 


represents a Hilbert-Schmidt class perturbation of the matrix J_|_ with the constant 

coefficients a(n) = 1 and 6(n) = 0. In this case we consider J_|_ as an operator acting 
in the standard space of one-sided sequences In its turn, the (sp)-condition means 
that the related spectral measure da has an absolutely continuous component supported 
on [—2,2]. Moreover, the spectral density a'{x) with respect to the Lebesgue measure 
satisfies an explicitly given integral condition, which in particular means that a'{x) ^ 0 
a.e. on this interval. Besides that, the measure may have at most countably many 
mass points (the set X) outside of the given interval. Again, the corresponding set X 
satisfies an explicitly given condition, which in particular means that the only possible 
accumulation points of this set are the endpoints ±2. Finally, note that there is no 
restriction on the singular component of the measure da on the interval [—2,2]. 

Later, the authors jointly with David Damanik generalized their result on the case 
of perturbations of periodic Jacobi matrices. To state this theorem we need a couple of 
definitions. 

We define a distance between two one-sided sequences 6 = {6(n)}n>o and 6 = 
{6(n)}n>o from by 

dist^(6,6) = dist^(6, 6) := ^ |6(n) — 6(n)|^??^”, ?7G(0, 1). (1.2) 

n>0 
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The distance dist(J+, J+) between two Jacobi matrices is defined via the distances 
between the generating coefficient sequences. 

Let J{E) be the isospectral set of periodic two-sided Jacobi matrices with a given 
spectral set Li C M. The distance between J_|_ and J{E) is defined in a standard way 

dist(J+, J{E)) = inf{dist(J+, J+) : J G J{E)}, 

O O 

where J+ is the restriction of a two-sided matrix J on the positive half-axis. 

(4) Damanik-Killip-Simon Theorem (DKST) [9]. Assume that J+ is a Jacobi matrix 
and let da be the associated spectral measure. The following are equivalent: 

(opp) Let S+ denote the shift operator in l'^. Then 

dist^{{Sir J+Sl, J{E)) < oo. (1.3) 

n>0 


(spp) The measure da is supported on Li U X, and moreover 



(1.4) 


Remark 1.2. Note that (1.3) means that the shifts of the given Jacobi matrix 

O 

approach to the isospectral set J{E), but possibly not to a specific element J of this 
set. In the same time (1.4) looks as a straightforward counterpart of condition (1.1). 

Remark 1.3. Let us point out that the spectral set of any periodic two-sided Jacobi 

O 

matrix J is a system of intervals of a very special nature: the system of intervals 
E = [bo,ao] \ U|^;^(aj,b j) represents the spectrum of a periodic Jacobi matrix if and 
only if Li = T“^([—2, 2]), where Tm{z) is a polynomial with only real critical points, that 
is, Tr{c) = 0 for c G R, and its critical values Tjn{c) obey the conditions \Tm{c)\ > 2. 
Actually, the Damanik-Killip-Simon Theorem was proved under an additional regularity 
condition \Tm{c)\ > 2 for all critical points c. In this case the degree m = g + 1- 

The paper [7] reviews recent progress in the understanding of the class of so-called 
finite gap Jacobi matrices and their perturbations. In the end of the article the authors 
posed the following question: “Is there an extension of the Damanik-Killip-Simon the¬ 
orem to the general finite system of intervals E case?” In the present paper we solve 
completely this problem, see Theorem 1.5 below. Naturally, this question was posed as 
soon as the original Killip-Simon theorem was published or even presented or proved. 
From this point of view [7] is just an explicit recent reference. 

Finite gap Jacobi matrices were discovered in the context of approximation theory 
[2, 3], [5, Chapter X]. They became especially famous because of their relation with the 
theory of integrable systems, for historical comments we would refer to [24] with many 
references therein. But the true meaning of this class was significantly clarified recently 
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by C. Remling: for a system of intervals E the finite gap class J{E) consists of all limit 
points of Jacobi matrices with an essential spectrum on E, having this E as the support 
of their a.c. spectrum. 

(5) Remling Theorem [32]. Let E he a system of intervals. Let J+ be a Jacobi matrix 
with the generating coefficient sequences {a(n),6(n)} such that its spectrum cj(J+) = 
E Li X, where X is a set of points, which accumulate only to the endpoints of the 
intervals, and ct'(x) / 0 for a.e. x G E. 

a{n) = lim a(n + mfc), b{n) = lim 6(n + mfc), 

m/g—>-H-cx) m/e^+OD 


for all n G Z, then the corresponding two-sided Jacobi matrix J belongs to J{E). 

Note that the system of shifts J+S'!|!:}„>o forms a precompact set in the 

compact-open topology (generated by the distance (1.2)). 

For E = [bo, ao]\u|^^(aj, bj) the class J{E) represents a ^r-dimensional torus, which 
can be parametrized explicitly. 

(6) Baker-Akhiezer parametrization for the elass J{E), see e.g. [40, Theorem 9.4]. For 
a € M^/Z^ let 


_20{a +i^ + a)e{a-i^ + a) - e{a- p +a) 

A[a) = a - ^ ^ -, B{a) = h + d^\a- 


9{a J- aY 


9{a J- a) 


where 

9{z) = 9{z, ^)=Y1 z G 

neZ9 

with the following system of parameters depending on E: 


(1.5) 


• U is a symmetric g x g matrix with a positive imaginary part, Imll > 0; 


• a G is an appropriate shift; 

• /i G M^/Z^ and ^ G are certain fixed directions of discrete and continuous 
translations on the torus M^/Z^, respectively; 

• a > 0 and 6 G M are normalization constants. 

Then J G J{E) if and only if 

a{nY = A{a — fin), b{n) = B{a — fin), (1-6) 

O 

for some a G M^/Z^. In this case we write J = J(a). Thus, 

J{E) = {J{a) : a G (1.7) 


Definition 1.4. For an arbitrary finite system of intervals E, we say that a Jacobi 
matrix J_|_ belongs to the Killip-Simon class KS(i7) if for some X the corresponding 
spectral measure da is supported on U and obeys (1.4). 
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Theorem 1.5. J+ belongs to KS(£') if and only if there exist eo(n) G and 

^ain) G /+, €b{n) G l\ such that (ef. (1.6)j 

n n 

a{nf = A(^ea{k) - lan) + ea{n), b{n) = B(^ea{k) - fxn) + eb{n), (1.8) 
k=0 k=0 

where A{a) and 13(a) are defined in (1.5). 

Remark 1.6. In the one interval case the functions A and B are constants, e.g. if 
E = [—2,2], then A = 1 and B = 0 and we obtain the original Killip-Simon Theorem. 

Remark 1.7. It is easy to see that a Jacobi matrix of the form (1.8) satisfies (1.3), see 
Lemma 7.2. Moreover, from our explicit formulas one can give immediately a suitable 
approximant for (S’f)^J+Sf, this is J(an) G J(E), an = Y^^=o^a(k) — frn; or conclude 
that, if the series fd := conditionally converges, then the coefficients of J+ 

approach, in fact, to the coefficients of the fixed element J(/3) G J(E), 

a^{n) — A(f3 — jan) —0 and b(n) — B(I3 — fin) 0, where n —oo. 

The representation (1.8) contains a certain ambiguity, for the reason see Remark 7.1. 


1.2 Basic ideas of the method and the structure of the paper 

The proof of DKST was based on two things: 

(i) Magic formula for periodic Jacobi matrices 


(ii) Matrix version of the Killip-Simon theorem 

The first one is the following identity. Let S be the shift in the space of two sided 
sequences P. E = [bo,ao] \ Uj^;^(aj, bj) = T~^^([—2,2]), then 


Tg+i(J) = 5^+1 + 


O 

for all J G J{E). The last matrix can be understood as the (s' -|- 1) x (g + l)-block 

O O 

Jacobi matrix with constant block coefficients A{n) = Ig+i and B(n) = O^+i. 

Now, for J+ the matrix Tg+i(J+) is a (2g + 3)-diagonal matrix, or, also a one-sided 
(s J- 1) X (s + 1) Jacobi block-matrix, see survey [10], 


Tg+l{J+) 


B(0) A(l) 

A(l) B{1) A(2) 


Such matrix has a spectral (s J- 1) x (s + 1) matrix-measure, say dE. According to [9] 
the matrix analog of (1.1) is of the form 

logdetH'(y)|V4^^dy+ ^ yi/fc-4 < oo, (1.9) 

yk& 
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as before [—2, 2] U Y is the support of cE. On the one hand this condition can be 
rewritten by means of the spectral measure da of the initial Jacobi matrix into 
the form (1.4), y = Tg^i(x). On the other hand, due to the matrix version of the 
Killip-Simon theorem, (1.9) is equivalent to Tg+i{J+) — + (5!()^+^) belongs to the 

Hilbert-Schmidt class. This is a certain bunch of conditions on the coefficients of J+, 
but we should recognize that extracting from this simple-looking condition the final one 
(1.3), is a very non-trivial task. 

Our first basic observation is the following. 

Lemma 1.8. For a system of intervals E there exists a unique rational funetion A.{z), 
A(oo) = oo, such that 


E = [bo,ao] \ lJ(aj,b,') = A-1([-2,2]), 
i=i 


and Im A( 2 ;) > 0 for Imz > 0. 


Proof. Let 4^(2:) be the Ahlfors function in the domain C\ E. Among all analytic 
functions in this domain, which vanish at infinity and are bounded by one in absolute 
value, this function has the biggest possible value Capa{E) = \z'^{z)\z=oo (the so-called 
analytic capacity) [1]. As it is well known [30] 


1-T(2) 
1 4 '( 2 ;) 



z-hj' 


Then 


A(2) 


1 


4'(2;) 


^qz -|- cq -|- 


E 


i=i 


A. 

Cj-Z' 


where Xj > 0, j > 0, and 'l'(cj) = 0, Cj G (aj,bj), j > 1. 


( 1 . 10 ) 


( 1 . 11 ) 

□ 


Note that in this proof we represented ^{z) as a superposition of a function T : 
C \ LI —)■ B with the Zhukovskii map. Essentially, (1.11) is our generalized magic 
formula, though it holds of course not for Jacobi matrices. 


Jacobi, CMV and SMP matrices. Jacobi matrices are probably the oldest object in the 
spectral theory of self-adjoint operators. They are generated by the moment problem 


[4] 



( 1 . 12 ) 


In this problem we are looking for a measure da supported on the real axis, which 
provides the representation (1.12) for the given moments {sfc}fc>o. In this sense CMV 
matrices are related to the trigonometric moment problem, which corresponds to the 
same question with respect to a measure supported on the unit circle. Note that this 
problem is also classical [4], but corresponding CMV matrices are a comparably fresh 
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object in the spectral theory [33, 34], The strong moment problem corresponds to 
measures on the real axis in the case that the moments are given for all integers k. 
An extensive bibliography of works on the strong moment problem can be found in the 
survey [19], concerning its matrix generalization see [37, 38]. 

As usual the solution of the problem deals with the orthogonalization of the gener¬ 
ating system of functions, that is, the system 

-1 (- 1 )^ 2 
X 

in the given case. The multiplication operator by the independent variable in 
with respect to the related orthonormal basis we call SMP matrix (this is exactly the 
way of the appearance of Jacobi and CMV matrices in connection with the power 
and trigonometric moment problem, respectively). In another terminology they are 
called Laurent-Jacobi matrices [6, 12, 18]. Very similar to the CMV-case, this is a five- 
diagonal matrix of a special structure, say A+ = A^{da). We assume that the measure 
is compactly supported and the origin does not belong to the support of this measure. 
In this case our A_|_ is bounded, moreover is also a bounded operator of a similar 
five-diagonal structure (just shifted by one element!) 

Note that, by a linear change of variable, we can always normalize an arbitrary two 
intervals system to the form ci = 0, see (1.11), that is, 

= [bo,ao]\(ai,bi) = A-i([-2,2]), A(z) = Aq + cq - (1.13) 

Without going into detail, dealing with the structure of SMP matrices, we can formulate 
our second basic observation. 

Proposition 1.9. [13] Let A{E) be the set of all two sided SMP matrices of period two 

O 

with their spectrum on E (1.13). Then A G A{E) if and only if 

A(l) = Aol + Co - Ai(A)-i = 52 + 5 - 2 . (1.14) 

O O 

Remark 1.10. It is highly important in (1.14) to be hold that both A and (A)“^ are 
five-diagonal matrices. 

Naturally, (1.13)-(1.14) have to be generalized to the multi-interval case. This leads 
to the concept of GMP matrices (G for generalized), see the next subsection. However, 
even after such a generalization the result on spectral properties of (“some”) GMP 
matrices of Killip-Simon class would be interesting probably only to a small circle of 
specialists, working with the strong moment problem. The point is that GMP matrices 
are used here as a certain intermediate (but very important) object. In a sense, this 
is the best possible choice of a system of coordinates. We can try to clarify the last 
sentence. The standard point of view on J{E) is to associate it with the hyperelliptic 
Riemann surface = {{z,w) : = nj=o(^ “ ^j){^ ~ bj)}- Then J{E) corresponds 
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to the “real part” of the Jacobian variety Jac(9J_E:) of this surface see e.g. [24, 25]. 
Periodic GMP matrices, satisfying 

A(l) = 5®+^ + (1.15) 

for A( 2 :) given in (1.11), are most likely the best possible choice for a coordinate system 
on the affine part of Jac(9lE), at least in the application to spectral theory. 

Thus, the point is to go back to Jacobi matrices. Let da be compactly supported 
and 0 does not belong to its support. We can define the map 

F+ : SMP ^ Jacobi 

just setting J+((t) in correspondence with the given A^{a). If so, we can define (in a 
naive way) a discrete dynamical system {Jacobi flow on SMP matrices) by the map J7+, 
which corresponds to the following commutative diagram: 


SMP ^ 

SMP 


.^+i 

.^+i 

(1.16) 

Jacobi 

Jacobi 



where <S+J+ = 5!^ J+S'+. 

The third basic observation deals with the idea of getting properties of the class 
KS{E) from the corresponding properties of the class of SMP (or, generally, GMP) 
matrices using the above introduced dynamical system A^{n) = J7’°"'(yl+). 

In the next subsection we give formal definitions for GMP matrices and the Jacobi 
ffow on them, but probably we can already outline the structure of the current paper: 

Section 2. Probably this is a bit unusual, but we start with inverse spectral theory for 
periodic GMP matrices. We recall the functional model for finite gap Jacobi matrices. 
In this model each operator is marked by a Hardy space H‘^{a) of character-automorphic 
functions in the domain C\E, where a is a character of the fundamental group of this 
domain (2.1), so, as before, a G cf. (1.7). Here J{a) is the multiplication operator 

by the independent variable with respect to the basis {e"}nez (2.3), and {e“}n>o is an 
intrinsic basis in H'^{a). The point is that in this domain the inner function 'I'(z) and 
the fixed ordering C = {ci, ■.■,Cg} of its zeros generate another natural basis {fn}n>o 
in H‘^{a) (2.7). Thus, we obtain a new family of operators (the same multiplication 
operator in the new bases) 

A{E, C) = {A(a, C) : aG 

This is the collection of all periodic GMP matrices associated with the given spectral 
set E and a fixed ordering C of zeros of the Ahlfors function '^(z). The fact that ^'(z) 
is single valued (the character corresponding to this function is trivial) is responsible 
for the periodicity of an arbitrary A{a, C). 



Another characteristic feature of '^( 2 :) is its certain conformal invariance. Indeed, if 
w = Wj = then 'l'j(tc) := 'I'( 2 ;) is the Ahlfors function in the tc-plane. The given 
ordering C generates the specific ordering 


1 





1 

0 ,- 

Cl Cj 


1 

Cj_i — Cj 


and the multiplication by w is again a periodic GMP matrix (up to an appropriate 
shift). That is, 

S-^{cj - A{a, C))-^S^ € A{Ej, Cj), (1.17) 

where Ej = {y = ■ x G E}. Note that 0 = w(oo). Let us point out that 

the spectral condition (1.4) possesses the same conformal invariance property. Thus, 
passing from the e-basis to the f-basis in we payed a certain prize: J{a) is three 

diagonal and A(a, C) is a {2g + 3)-diagonal matrix. In the same time we essentially 
win, since (cj — has infinitely many non-trivial diagonals, but due to (1.17) all 

matrices (cj — A{a^ C))“^ are still {2g + 3) diagonal. For them (1.11) (in the chosen 
basis) is nothing but the magic formula (1.15). 

The Jacobi flow on A{E^ C) can be defined in a very natural way. Since S~^J{a)S = 
J(a — ^) is a shift by a fixed character g, we set JA{a, C) = A{a — p,C). As we see, 
this is just one, probably new, object in the family of integrable systems. 

As a result, thanks to this section we are well prepared to understand and describe 
the structure of GMP matrices, A G GMP(C), and the Jacobi flow on them, A(n) = 
J'°'^A, in the general case. This is done in the Sections 3 and 4, respectively. 

In Section 5 we work with the Killip-Simon spectral condition for two-sided Jacobi 
and GMP matrices. Let us explain this passage to two-sided matrices. Our dehnition 
(1.16) is naive for the following reason. In the transformation J+(n) = the 

O 

eigenvalues in the gaps start to move. E.g., in a generic case for an initial J+, which 

O 

corresponds to one of our fundamental operators J G J{E), the eigenvalues will cover 
densely the spectral gaps (aj,bj). Thus, corresponding to such measures A+(n) just 
can not be properly defined. The easiest way to explain that nevertheless our program 
is doable is the following: use two-sided Jacobi matrices and enjoy unitarity of the shift 
S in l‘^l^ We show that an arbitrary one-sided Jacobi matrix J+, with its essential 

O O 

spectrum on E, can be extended by a Jacobi matrix J_ = P-JP-, J G J{E), such that 
for the resulting two-sided matrix J 


{cj — J) ^ exists for all cj. 


(1.18) 


see Lemma 5.1. We can improve the diagram (1.16), see 1.28, using the map J = EA 
on two-sided matrices, see Definition 1.14. In Proposition 5.5 we describe its image, 
which consists of Jacobi matrices with the property (1.18). 

^One can actually work with one-sided matrices but still use methods related to two dimensional 
cyclic subspaces. 
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Using the block-matrix version of the Killip-Simon theorem, it is a fairly simple task 
to write the necessary and sufficient condition for A € GMP(C) with the spectral data 
(1.4) in the form 

A(A) — -)- 5'®+^) is in the Hilbert-Schmidt class. (1-19) 

Note that the relation between corresponding spectral densities of V(A) and A has a 
quite elegant form (5.7). 

Section 6. Condition (1.19) is equivalent to 

II±(A} < oo (1.20) 

for the Killip-Simon functional of the problem, which is basically the /^-part of the trace 
of (A(yl) — (5'“^®"’“^^ -|- 5'^"’“^))^, for the precise expression see (6.2). In the spirit of our 
third basic observation, we compute the “derivative” of this functional in the direction 
of the Jacobi flow, that is, the value 

5jH+{A) ■.= H+{A)-H+{JA), 

see Lemma 6.1. This derivative represents a finite sum of squares! Now, we can rewrite 
(1.20) as the “integral” thus, get certain /^-properties. 

Note that they are already more related to the Jacobi matrix J = !FA, than to the 
given GMP matrix A itself. Nevertheless, all these conditions were given by means of 
the coefficients of A(^), not by the ones of A (or the system of iterates A{n), to be 
more precise). This is probably the hardest technical part of the work. To indicate 
the difficulty, we would like to mention the following. In [26] we found higher-order 
generalizations of Killip-Simon sum rules (relations between coefficients of J+ and the 
spectral measure da), for a single interval spectrum. But only for a very special family 
(related to Chebyshev polynomials of an arbitrary degree n), which was initially found 
in [22], we were able to convert the result of the form (1.19) to explicit relations on the 
coefficients of the given J_|_. Otherwise, each particular case becomes a reason for an 
interesting research, see e.g. [21, 16, 36]. Moreover, a nice looking general conjecture 
was recently disproved by M. Lukic [23]. By the way, for a highly interesting new 
development in this area see [15]. So, in this section we prove Theorem 1.20; practically, 
this is already a parametric representation for coefficients of Jacobi matrices of KS(K). 

Section 7. In this section we finalize the parametric representation for Killip-Simon 
Jacobi matrices associated to an arbitrary system of intervals E, that is, we prove the 
main Theorem 1.5. In the end of this section we demonstrate implicitly our last basic 
for this paper observation that the spectral theory in the spirit of [8] could be more 
powerful than the classical orthogonal polynomials approach [4, 35], see Subsection 7.2 
and especially Remark 7.3. Explicitly this was demonstrated in [28, 41, 27], as well as 
in Section 2 of the current paper. At the moment we are not able to present a theory 
of spaces of vector bundles, which corresponds as model spaces to Jacobi matrices of 
Killip-Simon class (in full generality) even in a finite gap case. 

Basic facts with respect to one sided GMP matrices are given in the appendix. 
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1.3 GMP matrices and Jacobi flow on them in solving the Killip- 
Simon problem 

In this subsection we give formal definitions for the named objects so that in the end 
of it we are able to state Theorem 1.20. This is the main ingredient in our proof of 
Theorem 1.5. 


• Let {sn} be the standard basis in Depending on the context, is the set of 

square-summable one-sided sequences or the subspace of spanned by {en}n>o- 
In the last case := and P± : ^ l\. are the orthogonal projectors. 

• Let {(5fc}|^Q denote the standard basis in the Euclidian space 

• By T* we denote the conjugated operator to an operator T, or the conjugated 
matrix if T is a matrix. In particular, for a vector-column p G 0+^, {p)* is a 
{g + I)-dimensional vector-row. Consequently, the scalar product in can be 
given in the following form (p, ^ = {q)*p. 

• The notation T~ denotes the upper triangular part of a matrix T {excluding the 
main diagonal), respectively T~^ := T — T~ is its lower triangular part {including 
the main diagonal). 


GMP matrices form a certain special subclass of real symmetric {2g -|- 3)-diagonal ma¬ 
trices, g > 1. First of all, the class depends on an ordered collection of distinct points 
C = {ci,...,Cg}. That is, if needed we will specify the notation GMP(C). We will 
define two-sided GMP matrices, but their restrictions on the positive half-axis will be 
highly important. 

Definition 1.11. We say that A is of the class A if it is a (s' -|- l)-block Jacobi matrix 


A*{p_i) B{p_i) A{po) 

A*{po) B{-po) A{pi) 


( 1 . 21 ) 


such that 


p = (p,^ G A{p) = dgp*, B{p) = {qp*) +{pq*)"^ + C, (1.22) 


and 



We call {pjljgz the generating coefficient sequences (for the given A). 


(1.23) 
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Definition 1.12. Let S be the shift operator Sen = Cn+i- A matrix A E A belongs to 
the GMP class if the matrices {cfc —are invertible, and moreover S~^{ck — A)~^S^ 
are also of the class A, see (1.21)-(1.23). To abbreviate we write A E GMP(C). 

Remark 1.13. As it follows from the definition ||5“^(cfc — A)“^5^|| < oo. These 
conditions can be written explicitly as a certain set of conditions on the coefficients 
{pjljez of the initial matrix A E A, see (3.7). That is, in fact, A E GMP(C) if and 
only if it is of the class A for a certain ordered collection {ci,... ,Cg} and (3.7) holds 
for the generating sequences. This can be regarded as a constructive definition of GMP 
matrices, see Theorem 3.3. 

Let J be a Jacobi matrix with coefficients {a(n),6(n)}: 

Jcn = a(n)en_i + 6(n)e„ + a(n + l)e„+i, a(n) > 0, n E Z. (1-24) 

The two-dimensional space spanned by e_i and cq forms a cyclic subspace for J. Also, 
J can be represented as a two-dimensional perturbation of the orthogonal sum with 
respect to the decomposition 0 

J= 0a(O)(eo(-,e_i) 0e_i(-,eo)). (1.25) 

We have a similar decomposition for A E GMP(C) 

[A- 0 1 1 

A= Q 0 ||po||(eo(-,e_i) 0e_i(-,eo)), e_i = e_i, eo := n^P+Ae-i. (1.26) 

Definition 1.14. For A E GMP(C) the Jacobi matrix J = J-'A is uniquely defined by 
the conditions 

’’±( 2 ^) := ((«^± - ^;)“^e^i±i,e 00 i) = ((A± - z)"^e^i±i,e^i±i), a(0) = ||po||- (1-27) 

2 2 2 2 

Definition 1.15. Let SJ := S~^JS. The Jacobi flow on GMP matrices is generated 
by the transformation 77, which makes the following diagram commutative 


GMP 


GMP 





(1.28) 

Jacobi 

—> 

Jacobi 



The corresponding discrete dynamical system (Jacobi flow) is of the form 

A(n0l) = 77A(n), A(0) = A. 

Essentially, it can be reduced to an open (input-output) dynamical system (4.12). 
The coefficients of the Jacobi matrix J = J^A are easily represented by means of the 
Jacobi flow acting on the initial A. Namely, 
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Corollary 1.16. Let J = FA and A{n) = J°'^A. In the above notations (1-24) 

a{n) = ||po(f^)||, &(n - 1) = (1.29) 

Now we can define the Killip-Simon class of GMP matrices. Let i? be a system of 
g + 1 disjoint intervals, E = [bo,ao] \ bj). Let A(z) = Ae{z) be the unique 

function, which was given in (1.11). 

O 

Proposition 1.17. A G GMP(C), generated by coefficients p = {p,q), belongs to the 
isospectral set of periodic matrices A{E, C) if and only if it obeys the magic formula 
(1.15). Moreover, the isospectral surface ZSe is given by, see (2.26), (2.27) and (3.2), 


Pg 


1 


Qg = -Co - Ao '^Pjqj, Afc(p) 
i=i 


Afc) h 1,g. 


(1.30) 


Definition 1.18. Let A G GMP(C). Let a± be the related spectral measures, that is, 


r±{z) 


da±{x) 
X — z ’ 


where r±(z) are given in (1.27). We say that A belongs to the Killip-Simon class 
KSA(£', C) if the measures a± are supported on K U X±, and both satisfy (1.4). 

The following theorem is essentially a consequence of the matrix version of the 
Killip-Simon theorem. 

Theorem 1.19. A G GMP(C) belongs to the Killip-Simon class KSA(£', C) if and only 
if the differenee Ae{A) — belongs to the Hilbert-Sehmidt class. 

However, the next statement is already highly non-trivial. Practically, it gives a 
parametrization of the coefficients of Jacobi matrices of Killip-Simon class with its 
essential spectrum on E. 

Theorem 1.20. For A G GMP(C), let A{n -|- 1) = J A{n), H(0) = A. Let {Pj(ra)}jez 
be the forming A{n) eoeffieient sequenees. A belongs to KSA(K, C) if and only if 

-pf\n)}n>o G 1%, - qf\n)}n>o G /+, (1.31) 

{AoP®(n) - l}n>o G ll, {>^o{poin),qo{n)) -h co}n>o G /+, (1.32) 

{Afc(po(n)) - Xk}n>o G ll (1.33) 

hold for all j = 0,g — 1 and all k = . ,g (ef. (1.32)-(1.33) and (1.30 )). 

To summarize, in this paper solving the Killip-Simon problem 


• we introduce GMP matrices as possibly the best coordinate system on the Jaco- 
bians of hyperelliptic Riemann surfaces associated to finite band operators; 
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• we introduce and study the Jacobi flow on GMP matrices as one more important 
object in a rich family of integrable systems; 

• our study is based essentially on the Damanik-Killip-Simon theorem on Hilbert- 
Schmidt perturbations of Jacobi block-matrices with constant coefficients; 

• we follow the ideology of application of analytic vector bundles in spectral theory, 
explicitly in Section 2 and implicitly in Section 7. 


2 Functional models for J{E) and A{E^Q). Jacobi flow on 
periodic GMP matrices 

2.1 Hardy spaces and the class J{E) 

In what follows, we will use functional models for the class of reflectionless matrices 
J{E) in the form as considered in [39]. To this end, we need to recall certain special 
functions related to function theory in the common resolvent domain Q = C \ E for 
J G J{E). Note that in this case, E can be a set of an essentially more complicated 
structure [17, 29, 43], than a system of intervals. 

Let ED/r ~ C \ iE be a uniformization of the domain fl. It means that there exists 
a Fuchsian group F and a meromorphic function 307 = 3 for all 7 G F, 

such that 


Vz G C \ F; 3C G D: 3 (C) = 2 ; and 3(Ci) = 3 (C 2 ) ^ Cl = 7 (C 2 )- 

We assume that 3 meets the normalization 3 ( 0 ) = 00 , (C3)(0) > 0- 
Let F* be the group of characters of the discrete group F, 


F* = {a\ a : F —> R/Z such that 0 ( 7172 ) = 0 ( 71 ) + 0 ( 72 )} 


Since F is formed by g independent generators, say the group F* is equivalent 

to R5/Z5, 


a~{a(7i),...,a(7g)} GR^/Z^. (2.1) 

Definition 2.1. For o G F* we define the Hardy space of character automorphic func¬ 
tions as 


H\a) = Hl{a) = {f e : / o 7 = 7 g F}, 

where denotes the standard Hardy class in B. 


Fix zo G and let orb(Co) = = { 7 (Co)} 7 Gr- The Blaschke product hz^ with 

zeros at 3 ~^( 2 :o) is called the Green function of the group F (cf. [39]). It is related to 
the standard Green function G{z, zq) in the domain H by 
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The function bzQ is character automorphic, that is, bzQ 07 = e^'^^^^obzQ, where fizo ^ T*. 
For bzQ we fix the normalization b^Q(O) > 0 if zq 7 ^ 00 and ( 3 b)(0) > 0 for the Blaschke 
product b related to infinity. 

We define = k°‘{C,Co) as the reproducing kernel of the space that is, 

</,A:^o> = /(Co) 

Remark 2.2. Let us point out that in our case this reproducing kernels possess a 
representation by means of 9 functions associated with the given Riemann surface [14]. 
As already mentioned, /c" has sense in a much more general situation, say, domains of 
Widom type. Although, generally speaking, they can not be represented via 9 functions, 
they still play a role of special functions in the related problems. 

Let k°‘{C) = /cq (C)) b(C) = bj(o)(- 2 )) and ^ We have an evident decomposition 

H^{a) = {z-}®bH\a-^^), e“ =(2.2) 

This decomposition plays an essential role in the proof of the following theorem. 


Theorem 2.3. The system of functions 

CiC) = b”(C) 


^^a-nM(O) 


(i) forms an orthonormal basis in H‘^{a) for n G N and 


(2.3) 


(a) forms an orthonormal basis in L?‘{a) for n € Z, 


where 

L\a) = {fGL^: / o 7 = 7 € T}. 


Proof Item (i) we obtain by iterating (2.2). A proof for (ii) in a much more general 
case can be found in [39, Theorem Ej. □ 


The following theorem describes all elements of J{E) for a given finite-gap set E. 

Theorem 2.4. The multiplication operator by i in L^(a) with respect to the basis {e“} 
from Theorem 2.3 is the following Jacobi matrix J = J{a): 


je" = a(n;a)e"_i -h 6 (n;a)e" -h a(n -h l;a)e"+i. 


where 


and 


a{n;a) = A{a — np), .4,(a) = ( 3 b)( 0)4 


' A:“(0) 
A:“+^(0) 


b{n;a)=B{a — np), B{a) 


3b(o) [ (n^(o) 

b'(0) \ A:“(0) 


(fe"+^/(0) l , (3b)^(0) 

A:«+/^(0) J b'(0) 


This Jacobi matrix J{a) belongs to J{E). Thus, we have a map from T* to J{E). 
Moreover, this map is one-to-one. 
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Remark 2.5. Using the representation of the reproducing kernels via 6 functions, see 
Remark 2.2, one gets A{a) and B{a) in the form (1.5). 

Remark 2.6. The following important relation is an immediate consequence of the 
above functional model 


S ^J{a)S = J{a — fi), Sen ■= Cn+i- 


(2.4) 


In particular, J{a) is periodic if and only if Nfi = Or* for a certain integer N. 


2.2 Class A{E, C) and Jacobi flow 

Now we turn to the functional model for A{E,C). The rational function A{z) and 
the single valued function '^(z), z G C\ E, were defined in (l.lO)-(l.ll). Let us list 
characteristic properties of ^(z): 

(i) |T| < 1 in U and I'l'l = 1 on Li, 


(ii) 'L(oo) = 'L(cj) = 0, I < j < g, otherwise '^{z) A 0- 


All this implies that 

Therefore ^(^(C)) = b(C) 0^=1 t’cj(C)- In particular, g, + Ylj=i = Or*. 

Let us fix G D such that 3 (Ci) = Cj and 7 j(Cj) = Cj for the generator 7 ^ of the 
group T. In order to construct a functional model for operators from A{E,C), we start 
with the following counterpart of the orthogonal decomposition (2.2): 

H\a) = {kl ^H\a) = {f^ ® • • • ® {f“} ® (2.5) 


where 


^-7TiaiA)m 

j!Qf _ _j!Q: _ _ 1 (,2 

V,(Ci) 


VC "(C 2 ) 

Theorem 2.7. The system of functions 


)•••) Ig 


V^A:“+^‘( 0 ) 


( 2 . 6 ) 


fn = fn(C;ci,...,Cg) = n = {g + l)m + j, j €[0,...,g] (2.7) 

(i) forms an orthonormal basis in H‘^{a) for n G N and 

(ii) forms an orthonormal basis in L?‘{a) for n G Z. 


Proof. Item (i) follows from (2.5) and for (ii) we have to use the description of the 
orthogonal complement L?‘{a) QH‘^{a), see [39]. □ 
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Similarly as we had before, this allows us to parametrize all elements of A{E, C) for 
a given E by the characters of F*. 

Theorem 2.8. In the above notations the multiplication operator by 3 with respect to 
the basis {f“} is a GMP matrix A{a; C) G A{E, C). Moreover, this map F* ^ A{E-, C) 
is one-to-one up to the identification {pj,qj) {—Pj,—qj) in A{E;C) ,0<j<g-l. 

Proof. We claim that the structure of the matrix is fixed by the choice of the orthonormal 
basis. In particular we need to check that, under the normalization (2.6), Pj{oi) and 
qj{a) are real. 

For fij = a - Yli=i hck, we have 


pAa) = = (l.3)(0) n bc.(0)e-"ftG) . (2.8) 

Since A:^(C) = k^{C) for all fi G F*, we get 

k^iCj) = kW) = kHljiCj)) = 

Therefore, Cj) is real. Note that the square root of 

g-27ri/3(7j) -g ^}^g niultiplier ±1. 

To compute the entries of the matrix B{p), and actually to show its specific struc¬ 
ture, we use a standard trick related to reproducing kernels. Let 0 < m < n < g. 
Then 

g-i7r(/3„(7„)-/3,„(7„))^/3n 

(sG, C> = (jbc™ ■ ■ ■ - =—— , fi ). 

Denoting for a moment the first function in the scalar product by f, since f(Cm) = 0; we 
can continue with 


= (3f-(l’3)(0)f(0) 



(b3)(0)K0) 


(b3)(0)(bc^ • • • 

V^fcJ(C„) HCm)k^-^n(o)^k^-{Cm) 

where, as before, Pn{oi) is of the form ( 2 . 8 ) and 


Pn{o:)qm{a) 


qm{a) : = 


) yp+F( 0 ) 

(bc2 ... bc_J(0)b(U)fe/5-+^(0)^A:g:(C™) 


Similarly, we get the representation for the diagonal terms. 
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If a periodic A G A{E,C) is given, we introduce its resolvent function r^(z). As 
usual it can be expressed by means of the transfer matrix, see Theorem 2.13 below. 
After that we define a exactly as in the Jacobi matrix case, see e.g. [39, (2.3.2)-(2.3.3) 
and Theorem A]. □ 

Definition 2.9. We define the Jacobi flow on A(E; C) as the dynamical system gener¬ 
ated by the following map (see (2.4), (1.28)): 

JA{a) = A{a — n), a G T*. 

We can describe this operation in a very explicit form. 

Lemma 2.10. Let 

^ Fsin,^ cos 4' ] p g, 

COS (j) — sm (j) 

Let 0{a) be the unitary, periodic (^f + 1) x (^r J- l)-block diagonal matrix given by 

0{oi) ... j ... 0{ot), (2.10) 

where 


0(a) = ^ 

^ |_ 0 o{(p{a) 


, [sin((/)(a)) cos((/>(a))] = 


[pg-i{a) Pgja) 

Jpi_i{a) +pl{a) 


( 2 . 11 ) 


Then 


OA{a- C) := S-^0{a)*A{a- C)0{a)S = A(a + Cg, ci,..., c^.i). (2.12) 

Proof. In fact, we consider here the following elementary operation on the set of periodic 
GMP matrices; we switch the oder in the orthogonalization procedure of the family of 
reproducing kernels from ..., A:“} to {A:"^, • • •, A:", A:"^}- Let {f“} be the new 

orthonormal system in Lf^a). It is easy to observe that 

t’c,(C)fn(C)=C^(C;C3,Ci,...,C3_i). 

That is, up to the shift, we derived a GMP basis of the form (2.6), but with the new 
ordering {cg, ci,..., Cg_i) and the new character a + //cg- Note that passing from one 
to another basis in a two dimensional space is a rotation o{(f). Thus the matrices of 
multiplication by 3 with respect to {f()(C; ci,..., c^)} and {f^ (C; Cg,..., Cg_i)} are 

related by (2.12). It remains to compute the angle f by means of A(q;;C). Since 
A{a + ticg',Cg,... ,Cg_i) G A, we have 


0 0 
Pg-i{a) Pg{a) 



That is, pg-i{a) cos 4> — pg{a) sin f = 0. Since Pg{a) > 0 we obtain (2.11) for o{(f) of 
the form (2.9). □ 


18 



Theorem 2.11. In the above notations 


JA{a\ C) = 0°^A(a; C) (2.13) 

Proof. We use (2.12), having in mind that = —p and that after all permuta¬ 
tions we obtain the original ordering C. □ 

The next lemma allows us to estimate components of the vector f", j = in 

its decomposition with respect to the basis {e“}n>o- 

Lemma 2.12. Let f" = YlT=o^ki^)^k- '^hen 

\Fi{a)\<C{E)rj'^, j = 0,...,g, (2.14) 


where l>g> max{|b(Ci)|, • • •, |l’(Cg)|}- 

Proof. First of all, we note that CfE) < ||A:"^|| < C{E) uniformly on a G F*. Also 
|l’c„(Ci)l — dehnition (2.3), |e^(Cn)| < c{E)g^. Since 


(efc, n 


= / __ 

ni=\K 

nf-i f. 




(Cn) nU,l^n (Cn) ' 

^ feg+^-(0)e^(Cn) 


ns (0) Si (O)feS^^" (Cn) nCl^n (Cn) 


and this is e ||F'^(q;)^ we get (2.14). 


□ 


2.3 Transfer matrix 

In this subsection we discuss briefly the direct spectral problem of the class A{E,C). 
For a vector v = p,q, we use the following notations 

(ufcu)* = [uo ... Vg-k], (dkv)* = [vk ... Vg]. (2.15) 

Recall, is the standard basis in We define inductively upper triangular 

matrices Mj’s of dimension [g + 1 — j) x (g + 1 — j) such that 

=M{p):=Mo= ° Ei-pqg + qPgWg (2.16) 

and 

Mj= ° +{-Ujpqg-j+Ujqpg-j)5*g_j, j>l. (2.17) 
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Theorem 2 . 13 . Let 


Roo(^) 

Bog{z) 


'{{Bo-z) ^p,p} 

{{Bo-z) ^5g,P)] 

Tlgo{z) 

Bgg{z)_ 


{{Bo - z)-^p,6g) 

0 

1 


(2.18) 


and r^(z) = ||p|p((A+ — z) ^eo,eo). Then the shift by one block for a one-sided GMP 
matrix yl+ i-)- see ( 2 . 20 ), by means of the spectral function has the following form 


r+{z) = 


'^ii{z)r^l\z) +2ii2(^) 

2 I 21 (^;)r^^^ (^;) +2l22(^;) 


where 


2 l(z) := 


■2I11 

2112' 

(^) 

1 

0 

0 

1 

oto 

1 

0 

0 

.2I21 

2 t 22 _ 

Rog{z) 

Bgg 1 J 


(z). 


(2.19) 


Proof. We represent as a two dimensional perturbation of the block diagonal matrix 


= 


B{p) 0 
0 




+ 




and apply the resolvent perturbation formula. 

Note that in the dehnition (2.19) we use the normalization det 2 l( 2 ;) = 1. 
Definition 2 . 14 . Let p* = [p ( 7 ] G M^. The matrix function 

1 


( 2 . 20 ) 

□ 


0(2;, c; p) = I - -^pp*j = e o-.PP \ j = 


c — 2; 


0 -1 

1 0 


( 2 . 21 ) 


represents the so-called Blaschke-Potapov factor of the third kind with a real pole c 
[31]. Note that (pp*))^ = 0. A specific factor related to infinity we introduce in the 
form 

ro -p 


a{z-,p) = 0(2:, 00; p) = 


1 z—pq 
Lp p J 


( 2 . 22 ) 


Theorem 2 . 15 . Let p* = [pj qj]. The matrix function 21 ( 2 ;), given in (2.19), pos¬ 
sesses the following multiplicative representation 


21(2:) = a(2:,ci;po)a(2;,C2;pi)... 0(2;,c^;pg_i)a(2;;p^). 


(2.23) 


Proof. We use the representation (2.16) and definitions (2.18), (2.22) to get 21 ( 2 ;) = 
2 lo( 2 ;)a( 2 :; Pg), where 


2 lo(^) = /- 


{{Mi-z) ^uip.uip) {{Mi-z) ^uiq,ui^ 

((Ml - z)-^uip, ui^ ((Ml - z)-^uiq, ui^ 


(2.24) 


Recall, that Ujp,Ujq were defined in (2.15). Then, we use one after another (2.17) and 
definitions ( 2 . 21 ) to get 


2 tj—1(2:) — 2lj(2:)ci(2;, Cg,_|_i_j; p^_j). 
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where 


. X r \{{Mj - z) ^Ujp,Uj^ {{Mj - z) ^Ujq,Ujp)'\ 

^ ^ [{{Mj - z)-^Ujp,Uj^ {{Mj - z)-'^Ujq,Uj^\ 

That is, we obtain (2.23). □ 

Definition 2.16. Let A G A{E,C). Then the product (2.23) is called the transfer 

O 

matrix associated with the given A. 

The role of the transfer matrix is described in the following theorem. 

O 

Theorem 2.17. Let A G A{E,C) with the transfer matrix 21 ( 2 ;), given in (2.23), and 

O 

let A{z) := tr 21 ( 2 ;). Then the speetrum E of A is given by 

E = A-\[-2, 2]) = {x : A{x) G [-2, 2]}. (2.25) 

Moreover, A{z) is of the form (1.11), where 

9 

XoPg = 1 , Ao + Co = 0 , (2.26) 

j=0 

and Afc = Afc(p) := —Rescj,tr 21 ( 2 ;), i.e., 

k—2 g—l 

Xk = -tr {fl <^i^k,Cj+i;pj)pk-iPl-i} n <^i^k,Cj+i-,Pj)a{ck;pg)}. (2.27) 

j=0 j=k 

Proof. A proof of (2.25) is the same as in the case of periodic Jacobi matrices. The 
relations (2.26) and (2.27) follow immediately from (2.23). □ 

Proof of Proposition 1.17. First of all, we have a parametrization of A{E,C) by the 
characters T*. It is evident that, in the basis (2.7), multiplication by T is the shift 
5 's+i, Thus, the magic formula for GMP matrices corresponds to the 

definition (1.11). The relations (2.26), (2.27) imply the form of the isospectral surface, 

O O 

that is, (1.30). Conversely, if A(x4) = + 5"^+^, then A is periodic by Naiman’s 

Lemma [9, Lemma 3.4] (alternatively, see proof of Theorem 8.4). Therefore we can 
apply Theorem 2.17. □ 


Later, in Section 6, we will use another representation for qg. 

Lemma 2 . 18 . qg allows the following alternative representation 

9 k-2 3-1 r-Q Q- 

k=l j=0 j=k L 

Proof. From the second relation in (2.26) and (2.24) one has 


^3 + ^=0 = 


tr 2 to( 2 ;) Q j_ (i 2 ; = ^Rescfctr 21 o (2 


0 ^ 


which is (2.28). 
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3 GMP matrices, general case. 


We hope after Theorem 2.17, and especially (2.27), it would be easy to perceive the 
following notations. 


Notation 3.1. For k = 1,... ,g the following functions (polynomials) are given by 
A*, =Af(pT‘’. ■ ■ ■ .phTlpl^lp... ,p“) (3.1) 


k—2 g—l 

= -tr{n ci(cfc,Cm+i;p^^’^^^)PfcJ^/^(Pfc2i)*j JJ a(cfc,Cm+i;pW)a(cfc;pJ^)} 

m=0 m=k 

If, as before, pm^ = Pm"*"^^ = Pm for all m G [0,51] this notation is simplified to 

// 

Afc(p) = Afc(po,---,Pfc-i;Pfc-i,---,P 9) = -Rescfetr2l(z). (3.2) 


Lemma 3.2. Let A G A. Then the formal inverse to — A is well defined as soon 
as 0. Moreover, eaeh element of the inverse matrix is a rational function of 

the coefficients pj,pj+i of the two consecutive blocks of A with the denominator 
In partieular, the vector /(c^) such that (c^ — A)f{ck) = Ck-i, k = with the 

veetor-components 

fi = = {(/®(cs))„}».„ e C>+‘. 

obey fj{ck) = 0 for j 0 {-1,0,1}, 

(/(-i))o = ... = if^-%-2 = 0, {f^% = ... = (/«), = 0, 


and 


A?i,fc(/("'^)fc-i = 1, A#,(/(i))fc_i = 1, 

= (pj-^)*i nS=fc a(cfc, c,-+i; p(-^))pi-^^ 


k,...,g-l, 




(p 


'm ) 








m = 0,..., A: — 2. 


(3.3) 

(3.4) 

(3.5) 


Proof. We have a purely linear algebra problem. To find /(c^) we solve the system 


-A(p-i)/-i = 0, 

(cfc - R(p-i))/-i - Ao/o = 0, 

-A*(po)/-i + (cfc - R(po))/o - A(pi)/i = 4-1, 
-A*(Pi)/o + (cfc - R(pi))/i = 0, 

-A*{p2)fi = 0 . 


(3.6) 


It is worth to recall that B{p) is an upper triangular matrix M{p) up to a one¬ 
dimensional perturbation, see (2.16), and its main diagonal in this case is C, see defini¬ 
tion (1.22). For this reason all inverse matrices to (c^ — i?(p)) can be found exactly like 
in the previous section in terms of products of Blaschke-Potapov factors a{z, c; p). □ 
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Theorem 3.3. Let A G A. yl belongs to the GMP class if and only if the forming 
sequences {pj-,qj} satisfy the following conditions 




jez 


> 0, for allk = l,...,g. 


(3.7) 


Proof. The formal inverse to — A can be found explicitly, see the previous lemma. 
Moreover, solvability of the system (3.6) is equivalent to (3.3). Thus, from one side, if 
Cfc — A is invertible we have 

= l((cfc —A) < ||(cfc —A) ^||. 

That is, (3.7) holds. 

_I4. 

In the opposite direction, the formal inverse operators exists since ^ 0 for all k 
and j. The following estimation is very useful 


|a(cfc 




< , due to a(cfc, cu p) = e 


-pp ) 


(3.8) 


C2(9)\\Af 

It shows that every non-trivial entry ((c^ — A)~^ej,en) is bounded by ^ 

see (3.3)-(3.5). In this estimation Ci > 0 depends only on the infimum (3.7). Since the 
formal inverse has only [2g + 3) non-trivial diagonals, we proved that all — A are 
indeed invertible, i.e., ||(cfc — A)“^|| < oo. □ 


4 Jacobi flow, general case 

Let us mention once again that Theorem 2.11 gives already a certain hint for a con¬ 
structive definition of the Jacobi flow. It will be defined via the unitary transformation, 
which after g rotations and one shift, maps GMP(ci,..., Cg) into itself. The first rota¬ 
tion creates the matrix A, which belongs (up to a suitable shift) to GMP(cg, ci,..., Cg_i) 
class. Then we create a matrix of the class GMP(cg_i, Cg, ci,..., Cg- 2 ), and so on... On 
the last step (making the shift) we get the required Jacobi flow transform, see (4.5). 
Having in mind (2.10) and (2.11), we give the following definition. 

Definition 4.1. We define the map 

O : GMP(ci,C2,...,Cg) GMP(cg,ci,...,Cg_i) 

in the following way. Let O = Oa be the block-diagonal matrix 
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where Ok are the (s' + 1) x (51 + 1) orthogonal matrices, see (2.9), 


Ok 



Then 

OA := SO\AOaS-^. 

It is required, but easy to check the correctness of this definition. Note 
p-entries of A = OA we get 

p} =p}_iCos(/)_i, l<j<g-l; 


Also, 



Qg 4 

1 

Sc 

II 

1 

* 

r (-1) (-1) , (-1) (-1)1 

qg_fPg_{+Cg 


~(0) ~(0) , ^ 


Po 

% Po +C9J 

1 

[ Qg-i Pg Qg Pg \ 


0 ( 0 - 

L Qg-iPg ' Qg Pg 'J 

Thus, the g-entries have the form 

Qo^pf'’ = - sin 0-i\/(Pg-i)^ + Qfp'h = (if-i\l{P^g-iY + {pTY- 

Our next definition is a counterpart of (2.13). 

Definition 4.2. We define the Jacobi flow transform 


J : GMP(ci,C2 , ...,Cg) GMP(ci,C2 , ...,Cg) 
by 

Let us note that 


5-(9+i)0(^)59+i = 0(5-i®+byi5S+i). 

This has an important consequence. 

Corollary 4.3. 

0{J°^A) = J°^{OA). 

Proof. Due to (4.5) and (4.6) we get 

J{OA) = = 0°(9+i)(5-(9+i)yi59+i) = 0{JA). 


(4.1) 
that for 

(4.2) 

(4.3) 

i). 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

□ 
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Let us turn to explicit formulas for the given transform. First of all, we note that 


JA = S-^UXAUaS, 
where Ua is a (5 + 1) x (51 + l)-block diagonal matrix 


(4.8) 


Ua = 


U{p-i) 

Uim) 


Lemma 4.4. In the notation above, see also (2.15), we have 


U{p)5o = 

WpW 


Uip)6k = 


1 


114 - 1 ^ 1114^1 


0 

II4^P 

-Pk-idkp 


l<k<g. 


(4.9) 


Proof. It follows from a step by step representation of the block U (p) as the product of 
orthogonal matrices, see (4.5), 


U{p) 


1 

0 

1 


0 

_ 1 

!_ 0 o( 4 )_ 


1 

0 


sin (fk cos (fk 


[Pk-i \\dkP\\_ 

114-1^1 


□ 

Theorem 4.5. Let ^(1) = JA and let ^(1), <?i^^(l)} be generating coefficient se¬ 


quences of A{1). Then 



= ||m II 

Pk^ 


Wrj II 

\\dkPj\\\\dk+iPj\\ 




Proof. We get (4.10) and (4.11) from (4.8) 


* 


pi^\l) 


(4.10) 

Pg-l(l) 



II 4 II 

{BiPj+i)pj+i,Pj+i) 

(4.11) 

Pg^\\Pj+i\\ II 4 + 1 P 

by Lemma 4.4. 

□ 


Remark 4.6. In view of Theorem 4.5 the Jacobi flow on GMP matrices can be related 
to an open (input-output) dynamical system, see (4.13). Let us fix a block-position 
j = 0, but vary n in A(n -|- 1) = JA{n). Then the coefficients related to the next 
block j = 1 are involved only in (4.11) and in a very specific way. If we define the two 
dimensional input by 


a-(n) = |4i(n)||, 6“(n) 


(R(pi(n))pi(n),pi(n)) 

WpiinW 
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and consider the parameters {po{n),qo{n)} as the internal state of the system, then the 
open dynamical system is defined by 


-1 

0 

c 

Po{n + 1)* 


UKPoin)) 

O' 

B{po{n)) 5ga^^{n) 

U{po{n)) 0 

[po{n + 1) 

-B(po(n + 1))_ 


0 

1 

a“(n)(5; 6“(n) 

0 1 


(4.12) 


and a°“*(n + 1) = \\po{n + 1)||. Note also that the output {a°“*(n + are 

the Jacobi parameters of J = TA (cf. (1.29), (4.16)) and the input is related to 
xhat is, this system represents the GMP transform on Jacobi matrices. 


A{0) 


A{1) 


A{2) 



(4.13) 


Theorem 4.7. Let ^1(0) := A € GMP(C), A{n + 1) = JA{n), n € Z, and e_i := e_i. 
Define 

Gm = Ua(o)SUa(i)S ■ ■ ■UA{m)Se-i, m > 0, (4.14) 

em-i = S-^U-I_^y ■ m < 0. (4.15) 

These system of vectors form an orthonormal system in P, with respect to which the 
following three-term recurrence relation holds 

Acm-i = a{m - l)em -2 + Km - l)em-i + a{m)em, (4.16) 

where a(m) and b{m — 1) are given by 

a{m) = ||po(W')|| = \JP^q\'^Y + ■ ■ ■ +P^g\mfi, h{m—l) = p^~^\m)K'~^\m). (4.17) 

That is, A with respect to {cm} is a Jacobi matrix J, moreover the transformation (4.8) 
corresponds to its shift S~^JS. 
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Proof. This relation for m = 0 follows basically from the definition. We consider the 
term related to e_ 2 , which is the most nontrivial in this case. We have 

A = S-^U^l_^^A{-l)UAi-i)S. 

Therefore l)C^yl(-i)eo- Using the block structure of Ua we obtain 

P.SAe., = C/-i_^)P_A(-l)U^(_i)eo = C/-i_,)e_i||po(-l)||. (4.18) 

Having in mind (4.16) and e_i = e_i, we formally define 

a(—l)e _2 ;= S~^P-SAe-i, a(0)eo ;= P+yle_i. 

Then, due to (4.18), 

a(-l)e_2 = 

which proves both (4.15) and the first relation in (4.17) for m = —1. 

We can write a similar relation for ^1(1). Using definition (4.8), we rewrite such a 
relation by means of the original A. As the result, we obtain (4.16) for m = 1, and so 
on. Simultaneously, we proved Corollary 1.16. □ 

Remark 4.8. In fact {en}^=-oo is a basis in P, see Proposition 5.5. That is, 

J = TA = F*AF, Fen = Cn, 


and F is unitary. 

5 Spectral conditions 

5.1 Killip-Simon spectral conditions for one- and two-sided Jacobi ma¬ 
trices 

First of all we mention the following 

Lemma 5.1. Assume that J+ is a one-sided Jacobi matrix with essential spectrum on 

O O 

E. Then it can he extended by a matrix J_ = P-JP-, J G J{E), such that each Cj 
belongs to the resolvent set (domain) of the resulting matrix (1.25). 

Proof. Let 


R{z) :=£*{J - z)-^£ = j 


f dT _ 

r-(z) ^ 

a(0) 

J X — z 

a(0) 

r+{z)-^_ 


-1 


(5.1) 


where £ : ^ P such that £ 

entries of R{z) we have 

1 


c_ 

C+ 


= c_e_i + c+cq. In particular, for the diagonal 


R-i-i{z) r-{z) 


^ +a{0fr+{z), - ^ 


Ro,o{z) r+(z) 


+ a(0)V_(z). (5.2) 
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If r+(cj) is zero or infinity, we choose J such that r_(cj) is regular, that is, r_{cj) ^ 0, 
r_(cj) 7 ^ oo. And vice versa, if r^{cj) is regular, we set r_(cj) = 0. In both cases 
i?_i^_i(cj) 7 ^ oo and -Ro,o(cj) 7 ^ 00 . Therefore the whole matrix (cj — J) is invertible. □ 

The spectral Killip-Simon condition can be formulated either in terms of measures 
a±, see Definition 1.18, or by means of the matrix measure dS from (5.1). 


Lemma 5.2. The measures a± both satisfy the Killip-Simon condition if and only if 
the matrix measure dS is supported on E UY and obeys 



(5.3) 


Proof. We note two properties of an arbitrary function F{z), which is analytic in the 
upper half-plan and has positive imaginary part. If such a function has a meromorphic 
extension in an interval (aj,bj) C M then its zeros and poles interlay. Secondary, F(z) 
is of bounded characteristic in the upper half-plane, and therefore 


r |logT(a;)| 

Jr l + x‘^ 


dx < 00 . 


(5.4) 


From the first property we get that all poles of the first and second functions in (5.2) 
satisfy the Killip-Simon condition in R\i?. Applying this fact once again we obtain that 
poles of R-i-i and iio, 0 ) that is the set Y, satisfy this condition. Similar observations 
show the opposite directions. 

With respect to the a.c. part of the measure we have 


r_(x) ^ a(0) 

-1 

1 " 

0 

_1 


r_(x) ^ a(0) 

a(0) r+(x)“^ 


0 c’''+(x) 


a(0) r+(x)“^ 


Therefore 


det S'(x) 


r_ (x) 

—rf^{x) -|- a( 0 ) 2 r_(x) 


2 

a'_{x)a'^{x). 


Applying (5.4) to r-{z) and —r_,_^(z) -|- a(0)^r_(z), we obtain an equivalence of the 
conditions for det S'(x) and cr'j.{x). □ 


5.2 Scalar and block-matrix spectral Killip-Simon conditions 

Theorem 5.3. Let A G GMP(C). Its spectral measure satisfies (5.3) if and only if the 
block Jacobi matrix A (A) belongs to the Killip-Simon class. 

Essentially, it follows from the Lemma 5.4 given below. We prove the corresponding 
lemma for a scalar measure a, assuming that cr{cj) = 0. Note, even if we start with 
an initial one-sided matrix such that a^{cj) > 0 for some j, due to the Lemma 5.1, 
we always can get a two-sided J such that S(cj) = 0 . Also, it is more uniform to set 
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A(z) = To pass to our case, where ^{z) is of the form (1.11), it is enough 

J '-j 

to send one of this Cj to infinity by a suitable linear fractional transform. 

So, let da be a scalar measure with an essential support on E = A“^([—2,2]) such 
that a{cj) = 0. We define the matrix measure dH by 


/ 


d^jy) . 

y-z ' 


^^-^W*{x)da{x)W{x), 


(5.5) 


where 


W{x) = 



(5.6) 


In other words, dH is the matrix measure of the multiplication by A(x) in 
with respect to a suitable cyclic subspace. Note that one can normalize this measure 
by a triangular (constant) matrix L such that L* j dE{y)L = /, that is, to choose an 
appropriate orthonormal basis in the fixed cyclic subspace. 


Lemma 5.4. Let H'(y) he the density of the a.c. part of the measure dE on [—2,2] and 
a'{x) be the density of da, respectively. Then 


det E'{y) 


nA(x)=y 

nLi A. 


(5.7) 


Proof Let {xi,... ,Xg} = A ^{y), y £ [—2,2], Then 
='(y) = 

A{x)=y 



'a-'(xi) 


r 1 

1 1 


A'(xi) 


Cl-311 

Cg—Xl 

w* 

o-'(xg) 

W, W:= 

1 

1 


i A'(Xg)j 


LCl-Xg 

Cg—Xg _ 


As it is well known 


detW = (-1) 


9 ( 9 - 1 ) Y[k<ji^k Xj) ^j) 

- ^k) 


(5.8) 


On the other hand. 


Therefore, 


y - A(x) = y 


nix 



That is. 


-A'(xk) = y 


nk^jjxk-Xj) 

Ujixk-Cj) 


and 


, nAck-xj) 

Afc = yj ^—7 -r. 

Vik^ji^k Cj) 


A'(xfc) 


. Uk^jjXk - Xj) Uk^ji^k - Cj) 
^ UjiXk-Cj) UjiCk-Xj) 
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Thus 


(5.9) 




Y\k<Mk - Xjf{Ck - Cjf 
Wk,Mk-XjY 




Combining (5.8) and (5.9), we obtain (5.7). 


□ 


Next, we prove the following general statement. 

Proposition 5.5. Let a two-sided Jacobi matrix J he such that Cj ^ cr^J). Then, up 
to the identification {pm,qm) {—Pm,—qm), there exists a unique GMP matrix A 
related to the fixed ordering C such that J = FA = F*AF. In particular F ■. ^ is 

unitary. 

We need to define a counterpart of a cyclic subspace (2.5) in the general case. 
Assume that J = FA, see (1.25)-(1.26). Recall F : P ^ P is the isometry 


F Cm. — Cr 


where Cm were defined in (4.14)-(4.15). In particular, Fc-i = e_i and TP+ = P+F, 
Fcq = Co = We note that 

{h={A- Cl)/ : f^ll, if, eo) = 0} = {h e ll : {h, eo) = 0}. (5.10) 


Thus, T*eo can be described by means of an orthogonal complement in the following 
construction. 

Let c 0 ct( J) and, actually, it is not necessary, but let c be real. Having in mind the 
previous paragraph, we define 

ll^^:= {h = {J-c)f : fell, (/,eo) = 0}. (5.11) 

Recall that r^{z) = ((J+ — 2:)“^eo,eo). 

Lemma 5.6. Let Ac = l'^ Q l%c- This is a one dimensional space, i.e., Ac = {xc}- 
Moreover, we can choose 

Kc = {J — c)“^(e_ia(0) sin (/? + cq cos ip), (5-12) 


where 


tanv? = tan(^(c) = r+(c), 


(5.13) 


including = f i/r+(c) = oo, that is, c is a pole of this function. In this notations 


Xr. = 


4(c) 




1 + r+(c)2 

Moreover, the following two-sided estimation holds 

min{a(0)2,l} , max{a(0)2,l} 

-^ < 7^ (c) < —2^- 


c + 


disR^(c, (t(J)) 


(5.14) 


(5.15) 
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Proof. If r+(c) 7 ^ oo, we have Xc = (J+ — c)~^eocos(p. Otherwise Xc is collinear to the 
corresponding eigenvector of J+. These prove (5.12), (5.13). 

Further, we have 


|>^c|| =((^+-c) "*60,60) cos"* v? = 


r+ic) 


1 + r+(c)2' 

which proves (5.14). We use (5.1). Since J dT, = I, we have 


1 


(|c| + 

Using (5.1), we obtain 
1 


< R’{c) = 




< 


(x — c)2 dist^(c, cr(J)) 


or 


(|c| + ||J||)2 

i?(c)-2 


< R{c) 


(|c| + ||J||) 


2 < 


rLjc) 

r- (c)^ 

0 

' F_(c) 
r_ (c)^ 

0 


0 


0 

’■+(0". 


i?(c) < 


1 


< 


dist^(c, ct(J)) 
i?(c)-2 


dist"^(c, (t(J)) 


Comparing the values in the lower corner of these matrices, we get 


1 + agr+(c)^ 

i\c\ + \\j\\r 

Thus, (5.15) is also proved. 


< r',(c) < 


dist^(c, (t(J)) 


□ 


Proof of Proposition 5.5. Defining Xc by (5.12), we obtain F*6o = \\ ^ci- In partic¬ 

ular. 




p!)°\o) = (Ae-i,eo} := (Je-ujA^) = 


K, 


Cl I 


a(0) sm(p{ci) 

p'{ci) 


Generally, we consider the ordered system of vectors 

Xci , . . . , X'cg ) 60 • 


(5.16) 


This system is linearly independent. Otherwise, there exists a nontrivial vector = 
{CjYj^Q such that 


0 = {{J-z) ^(6_iaor+(2;)-)-6o),XciCoH- ^ eoCg) = {{J+- z) ^6o,Xc,CoH-^eo^a) 

-{» + ■■■+ --''-'“U ,-! + r+{zK, 


r+(z) cos(^(ci) — sin0(ci) ^ ^ ^ r+( 2 ;) cos (/>(cg) — sin0(cg) ^ 


Z — Cl 


z-c„ 


COS(/>(ci)^ ^ ^ COS(p{Cg] 


Z — Cl 


-6 + • • • + 


Z -C„ 


4-1 + ^9 r+{z)- 


sin(/)(ci)_^ ^ ^ siiKpjcg] 


Z — Cl 


-4 + • • • + 


Z-Cn 


4-1 
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That is, r+(z) is rational and the corresponding has only a finite number of mass- 
points, which contradicts to the original assumption that a{n) > 0 for all natural n. 

In the spectral representation the system (5.16) corresponds to the vector-functions 


1 

ao sin(^(ci) 

1 

ao sin0(cg) 


'0 

X — Cl 

cos (^(ci) 

1 

0 

cos ^(Cg) 

1 

1 


from L^s- 
subspace 


Jointly with its orthogonal complement they form (2g'-|-2)-dimensional cyclic 


1 

0 ^ 

1 

1 ... 1 

1 

7 

7 


1- 

7 

0 

o' 

_l_ . . . _l_ 

X-Cg 

_ Co,5-1 _ 

. Co, a J 


(5.17) 


of the operator multiplication by A(x) in this space. Being ordered and orthogonalized 
in an appropriate way, it generates a GMP basis in The operator of multiplication 
by the independent variable with respect to this basis forms A G GMP(C). Moreover, 
its spectral matrix measure is dS, that is, the spectral measure of the initial J. □ 


Proof of Theorem 5.3. Glearly, the eigenvalue spectral condition on A corresponds to 
the eigenvalue spectral condition for A (A) of the Killip-Simon class matrices with 
asymptotically constant matrix-block coefficients. 

We see now that the cyclic subspace (5.17) of A(^) for a two-sided GMP matrix 
A represents indeed a simple two dimensional counterpart of the system (5.6). Subse¬ 
quently the matrix measure (E 2 {y) of A(yl) is of the form (5.5), i.e.. 


/ 


<E2{y) 

y - z 


We have 


-- Wf{x)di:{x)W 2 {x), W 2 {x) 

A(x) — z 


h 

Cl —X 


I2 

Cg-X 


h 


^2iy) 


A{x)=y 


A'(x) 


W 2 (x) 


A'(a:i) 


A'(a;g). 


W 2 , W 2 :=W(g)l2. 


Using det W 2 = det^ W, by a two-dimensional counterpart of (5.7), we obtain an equiv¬ 
alence of the Killip-Simon a.c. spectral conditions on A and A(^). □ 
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6 Theorem 1.20: the first step in a parametrization of 
coefficient sequences of the Killip-Simon class 

6.1 “Derivative” in the Jacobi flow direction 

Let us make the block decomposition of A(yl) in (gf + 1) x {g + 1) blocks 


A(A) = 


’-1 


tt)_i Oo 


•■o 


Wo 


( 6 . 1 ) 


where is a self-adjoint matrix and is a lower triangular one, i.e., 



(k) 

^0,0 ■ 

_1 


1- 

o 

O 

1_ 

Wfc = 



•) 

: 0 


L%,o • 

• • ^9,9_ 


Vg Q . . . Vg^g^ 


Due to the previous subsection and general results on Jacobi block-matrices of Killip- 
Simon class [9], the spectral condition (1.4) is equivalent to the boundedness of the 
following KS-functional 

H+{A) = ( 6 - 2 ) 

j>o 

where 

/i(t)j,rDj,t)j+i) = itr (t)*t)j + tul -t 0j+it)*+i) - (g-t 1) - logdetOjOj+i. 

Lemma 6.1. Let 


1 


6jH+{A) = -{AiJA)e.i,AiJA)e.i) - 1 


Then 


H+{A) = H+{JA) + 5.jH+{A). (6.3) 

Proof. Comparing JA = S~^U\A{A)UaS and U\A{A)Ua, we note that 5jHj^{A) 
adds to A) exactly that terms, which were omitted (beause of the shift) in the 

trace-like expression (6.2) for Pj^Uf^A{A)U a- Further, since Ua is of a block diagonal 
form, we have the following identities between the blocks of U\A{A)Ua and A(yl) itself; 

tr U*{pj)t)*t)jU{pj) = tr tr U*{pj)to‘jU{pj) = tr tUj, 

tr U*{pj)\3j+it)*^^U{pj) = tr 

Also, all Vj and U*{pj-i)t>jU{pj) are triangular matrices, and we have 

9 

= detOj = detU* {pj-i)\:)jU{pj). 

1=0 

After that, we arrive to the conclusion that the right and left hand side in (6.3) coincide. 

□ 
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6.2 “Derivative” in the GMP direction 

Let © be unitary in such that & ■. ^ l\. We denote by /C© = l\ © 0/^. For 

0„ = Ua(o)SUa(i) ■ ■ ■ UA{n)S we have 

^©„59+l = ® ©nAC^g + l = /C 59 +I © S^^^lCg-(g + l)Q^gg + l. (6.4) 

The system eo,... ,en, © = forms a basis in /C©„, see (4.14), and bq,. .. ,eg is 

the standard basis in ICgg+i. By ej = &jS^~^^e-i = ®jeg, j = 0, ... ,n, we denote a 
similar orthonormal system in 5^'''^/C5-(g+i)0^5g+i. 

Remark 6.2. We will need both, the standard shift for GMP matrices Ai := 
as well as the shift of A in the Jacobi flow direction, i.e., ^(1) = A. For this reason we 
have to use quite complicated notations. Recall that denotes the {k,m) entry of 

the the block X>j{n) in the matrix A(yl(n)), where A{n) = J°'^A. Then has the 

same meaning with respect to the shifted GMP matrix yli(n) = J°'^Ai = {J°^A)i = 
A{n)i, i.e.: W = OfcrlWi- 

Lemma 6.3. Let us define the partial sum for the functional H^{A) by 

1 

H+,niA) := ^{-(A(^(m))e_i,A(^(m))e_i) - 1 -logo5"^)(m)t)®(m)} 

m=l 

and let, as before, 

h{'Oo, tuo, Oi) = ^tr {OgOo + Wg + OitJ);} - (g + l) - log det Oof 1 . 

Then 

H+,n{A) + /r(Oo(n),tUo(n),Oi(n)) = h(Oo(0), tUo(O), Oi(0)) + R+,„(^i). (6.5) 

Proof. Let Pjc^ be the orthogonal projector on /C©^ 59 +i. We compute the trace 
of the matrix 

l^©„S9+l 

with respect to the two decompositions (6.4). In the first case we get 

n n 

^{A{A)em,A{A)em) = ^ (A(A(m))e_i, A(A(m))e_i) (6.6) 

m=0 m=0 

and 

9 

'^{A{A)Qnej,A{A)Gnej) = tr {Oo(f^)*Oo(n) + tuo(n)^ + 0i(n)0i(n)*}. (6.7) 

1=0 
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With respect to the second orthonormal system we get 
9 

^(A(A)ej, A(A)ej) = tr {OoOq + Wq + OiOl;} (6.8) 

j=0 

and 

n n 

^(A(A)em, A(A)em) = ^ (A(Ai(m))e_i, A(yli(m))e_i). (6.9) 

m=0 m=0 

Thus the sum of the expressions in (6.6) and (6.7) is the same as in (6.8) and (6.9). 
Next, we claim that for an arbitrary n 

n n 

vSo’(O)... 4J(0) • Yl (n)... vf^{n). (6.10) 

m=l m=l 

Recall that v‘fg{m) can be regarded as Vgg^\m)i, see Remark 6.2. We prove this by 
induction. Since 

uiJ(O)... u® (0) = det Oo(0) = det U*{p-i)\)oU{po) = ■ uJJ (1)... ^^^^.^(l), 

we have the following recurrence relation 

( 1 ) ■ ^00 ( 1 ) • • • 4 ? ( 1 ) = ^00 ( 0 ) • • • 4 ? ( 0 ) 4 ? (!)• 

That is, 

4o (0) • • • 4? (0) • 4? (1)4? (2) =44^ (1) • 4o (1) • • • 4? (i)4? ( 2 ) 

=44^(i)44^(2)-4?(2)---4?(2)> 

and so on... Thus, (6.10) is proved. 

In a combination of (6.6), (6.7) and (6.8), (6.9) with (6.10) we obtain (6.5). 


Theorem 6.4. If A) < oo, then H^{A) < oo. If 

Hj^{A) < oo and liminf/i(do(n), tDo(n), t)i(n)) = 0, 

n^oo 

then the Killip-Simon functional is finite, moreover (A) = i7+(T). 

Proof. We get II{A) < oo, iterating (6.3), morovere we obtain 

H+{A) = H+{A) + lim H+{J°"^A). 

m—yoo 

If (6.11) holds we can pass to the limit in (6.5) 

H+iA) = /i(oo, tt)o, oi) + i7+(5-(^+^)yl5^?+i). 

Iterating this identity we get that 77+(yl) is finite, and, in fact, 

H+{A) = H+{A)+ lim 

m—>-oo 

Therefore, by (6.12), 

H+{A) = H+{A) + lim + lim H+{J°"^A). 

m—>-oo m^oo 

Since both limits are nonnegative, we obtain 77+(A) = 77+(yl). 


□ 


( 6 . 11 ) 


( 6 . 12 ) 


□ 
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6.3 Proof of Theorem 1.20 

Lemma 6.5. Let A € KSA(i5^, C) and A{n + 1) = JA{n), A(0) = A. Then (1.31) and 
the first relation in (1-32) are satisfied. 

First we prove the following sublemma. 

Lemma 6.6. Assume that for sequenees fin and fin there are sequences Tn and fn such 
that 

G ll, (6.13) 

that is, all entries of the above matrix form l\-sequences. Assume in addition that there 
is rj > 0 such that for all n we have a priori estimations 

cos fin > v, cos fin > ??, - > Tn > V, - > Tn > f]. (6.14) 

T] T] 

Then - e^^'-}n>o G 1%. 

Proof. Directly from (6.13) we have 

{cOsV’n - cosfin}n>0 G and {Tn cos fin - fnCOsfin}n >0 G 
Then (6.14) implies {r^ — fn}n>o £ Now, we have another two conditions 
{TnSinfin - smfin}n>0 G and [sin fin - fnAnfin}n>C G l\. 


Tn 0 
0 1 


o{fin) - o{fi„ 


1 0 
0 fn 


Therefore, 


sin fin - Tnfn sin fin - fn (sin fin - Tn sin fin) 

belongs to that is, {sinV'„(l — Tnfn)}n>o & Thus, (r^ — 
sequence, as well as (r^ — 1) sin fin. Finally, since 


1) sin fin forms a 


sin fin - sin fin = Tn sin fin - sin fin - {Tn - 1) sin fin, 


both {sin'0„ — sin'0„}„>o and {cosV’n — cosfin}n>o are /^-sequences. □ 

Proof of Lemma 6.5. The first relation (1.31) follows immediately from Lemma 6.1. 

Let A = 0{A), see (4.1). We use tilde for all entries related to A and A(T) (6.1), 
respectively. The entries of A{n) we denote by {p^l^\n), q^^{n)} and we use a similar 
notation for the entries of A{A{n)) and A{A{n)). Due to Definition 4.1, 


[ Xop'i'’{n) 


0 


o((/>W(n)) =o((/.J i)(n)) 
Applying Lemma 6.1 to the matrix A, we obtain 


Xopf^ (n) 0 

[ 42(") 4!o(«), 


(6.15) 


{AoP®(n) - l}n>o G ll, ^ ^ 
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Similarly for the entries related to A we have {Aoiig°^(n) — l}n>o £ {'^o^g('^)}n>o G 

1“^. Thus, we can apply Lemma 6.6 with respect to (6.15). We get {sin(n) — 

sin(/)®(n)} belongs to l\. That is, {p^Z\^(n) — p^gli{n)}n>o G l\- 

Using (4.2), (4.3), we get similar relations for all others j’s. Using (4.4), we prove 
the second part of (1.31). □ 

Proof of Theorem 1.20. Lemma 6.6 implies that {v^g^ig_i{n) — l)sin())g ^^(n) form an 
/^-sequence, or, equivalently, see (3.1), 

{(A^i_^(n) - \g)p^~}l{n)}n>o G li- ( 6 . 16 ) 

Since p^_\^(n) may approach to zero, it does not imply yet that g{n) — A^} belongs 
to /i. Let us show that 


{{Ati^g{n) - Ag)gJ_i^(n)}n>o G 1%. (6.17) 

Since inf„ ^(g^~^^(n))^ + (p^~\^(n))^^ > 0, both (6.16) and (6.17) give us (1.33) for 


m = g. 

To this end, we note that 


cos (pg (n) 


(6.18) 


Indeed, by definition of the Jacobi flow 


U{p- 2 {n)) 


-^ 9,9 

* V, 


(- 1 ) 

0,0 

* 

* 


* V 


(- 1 ) 


(n + 1) = V-i{n)U{p-i{n)) 


the second from below entry in the last column in this matrix identity means exactly 
(6.18). Therefore, by Lemma 6.6, we get 


{^\g{n + 1) - A*^^g{n)}n>o G 


(6.19) 


Now, by (4.10) 

i^ti^gin) - Xg)pl~_^^{n) = -{At^g{n) - Ag)g^li^(n + (n)\\dg-ip-i{r 


7_i(n) 


In combination with (6.19) we have (6.17), and therefore (1.33) for k = g. 

The same arguments with respect to O^A, k = l,...,g — 1, in a combination with 
(4.7), give (1.33) for all other k. 
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Now we start with the conditions (1.31)-(1.33). Let us show that they imply the 
second relation in (6.11). Due to compactness we can choose convergent subsequences 

O O 

Po{nk)^P and qo{nk) ^ q. (6.20) 

o o o 

Let A be the periodic GMP matrix generated by {p,q). Passing to the limit in (1.32) 
and (1.33) along the subsequence {n^} we get the isospectral conditions (1.30), that is, 

O o o o o 

A G A(E,C). By the magic formula, 0 = /, to = 0, where 0,tD are blocks of A(yl). 
Recall that blocks of A(A) are formed by the coefficients of two consecutive blocks of 
A, see Lemma 3.2. Therefore, due to (6.20) and (1.31), we obtain 


O O 

lim tt)o(rafc) = to, lim Oo(nfc) = hm Oi(nfc) = 0. 

k^oo k^oo k^oo 


Thus, Iimfc^ooh(0o(nfc),tt)o(nfc),0i(nfc)) = 0. 

To show the first relation in (6.11), we evaluate the entries of A(yl)e_i. Let h{ck) = 
(cfc — y4)“^e_i. In notations of Lemma 3.2 we have 


{h^ ^\ck))i = 




{h^ = 


4-1 (pi-V)*i 


i-i 

JJo(cfc,Cj+i;p4^^)p!”^^ 


j=k 


UO) 


foT k < I and {h^ ^^(c;+i)); = -, where 


TT 


(0) ■ 
fc-1 
'(0) 
Pk-l 


k-2 


i=o 


4-V Pk-i 


<?-! 


j=k 

( 6 . 21 ) 

We note that due to the uniform estimations (3.7) and (3.8), conditions (1.31) and 
(1.33) imply that 

f A. 




- 1 GlX. 

J n>0 


( 6 . 22 ) 


For I < g, by definition we have 

i 

= ^opl~^\ri)q^i~^\n) + Xi+i{h^~^\ci+i))i{n) + Xk{h^~^'> {ck))i{n) 


k=l 


We substitute in this expression p^^li from (6.21). Then (6.22) and the first relation 
in (1.32) imply that differs from the sequence {wig^\n)}, see below, by a 

/^-sequence; here 


w 


(- 1 ) 

1,9 


(n) 


:=pj 


n — 


P 

j=0 


( 0 ) 

b 


n))p4(?^)+ (6.23) 
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I k-2 


o(cfc,Cj+i;p®(n))p®^(n) 

k=l jr=0 


(pLV(n))*j 

C/+1 - Cfc 


l-l 

JJ a{ck,Cj+i;pY\n))pY\n) 


j=k 


Using (1.31) once again, we can substitute in the last expression all p^^\n) by pj ^^(ra). 
After that, we note that the product 


k-2 l-l 

a{ck,Cj+i-pY\n))p[z^^{n){p[~_^^in))*] JJ o(cfc, c^+i; pj"^^(n)) 


j=0 


j=k 


is the residue of the matrix function 0 ^= 0 p^- ^\n)) at c^. It remains to use 
the identity 


i 


E 


Res cfc 0^=0 a{z, Cj+i ; pj- (n)) 

Z-Ck 


l-l 

Yla{z,Cj+i-,pY\n)) - I 

j=0 


(6.24) 


evaluated at z = c;_|_i and we get cancellation of the first and second line in (6.23), i.e., 
{wY\n)} G As result, we obtain {it';~^^(u)} G for all 0 < / < 5 ^. 

The diagonal entry require a little bit special consideration. To prove that (n)} 

belongs to we use 

g (-1) ~(o) 

(A(A)e_i,e_i) = +co - ^ Y) ’ 

k=i Pg ^-i,k 

the second relation in (1.32), and, instead of (6.24), a more involved identity (2.28) 
shown in Lemma 2.18. Similarly, one can prove that {vl^g\n) — 1} and {u®(n)}, for 
0 < / < S', form sequences. Thus, by Theorem 6.4, H+{A) = iL+(A) < 00 and, due 
to the matrix version of the Killip-Simon theorem A G KSA(ili, C). 

□ 


7 Proof of the main Theorem 1.5 

7.1 Prom GMP to Jacobi 

Assume that A G GMP(C). Let A(n) = J°'^A. Recall that the coefficients of 
the Jacobi matrix J = TA are given by (1.29) and P properties of the coefficients 
{p±i{n),q±i{n)^po{n),qo{n)} are given in Theorem 1.20. We consider the isospectral 
surface XSe given by (1.30), see (2.27), with the identification {pj,qj) = i—Pj,—Qj), 
j = 0,— 1. Note that this is a (/ dimensional torus, which we can parametrize by 
a G M^/Z^ according to Theorem 2.8. Moreover, by statement (d) of Moser-Uhlenbeck 
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Theorem [25, Theorem 4.7], for the given manifold 


sup ||(T(p)*r(p)) ^||<oo, T(p) 

{p}eX5B 


Note that evidently 


■ dM 

dpo 

QAi 

dqo 


dAi 

dpg-i 

dAi 


dAg 

dpo 

dqo 


dAg 

%' 

J 


0 < 


inf 


||(T(p)*r(p))-i||-i 


< sup ||T(p)*T(p)|| < oo. 

{p}€X(S£; 


(7.1) 


O 

We define a periodic GMP matrix A{an) generated by {p(an)} G ISe such that 

O 

dist(po(n),X<S£;) = dist(po(n), p(an)). (7.2) 


Using the standard Lagrange multipliers method, we can estimate the distance from 
Po(n) to the isospectral set in terms of sup ||(T(p)*T(p))“^||. Then, by (1.32), (1.33) 
and (7.1), we have 

A 

^dist^(po(n),p(an)) < oo (7.3) 

n=0 

and also, see (1.29), 

a(n)^ — A{an) G , b{n) — B{an) G 1“^- (7.4) 

On the other hand, by (1.31)-(1.33) and the uniform smoothness of the Jacobi flow 
transform (4.10), (4.11) 

O O 

dist(po(n + l),p(an - m)) < C{E, J){dist(po(n), p(an)) + dist(po(n), pi(n))}. 


That is, 

o o o 

dist(p(an+i),p(an - 1^)) < C{E, J){dist(po(n), p(an)) + dist(po(n), pi(n))} 

O 

+ dist(po(n + l),p(an+i))- 


Since 

jja - /3|| < Ci(L;)dist(p(a),p(/3)), 

(7.3) and (1.31) imply 

OO 

^ ||ea(n)||^ < oo, where eain) := an+i - {an - fi). 

n=0 

In combination with (7.4), we obtain (1.8). 

Remark 7.1. Of course in this proof it is not necessary to choose an as the best ap¬ 
proximation to po(u.), see (7.2). It is enough to have this distance under an appropriate 
control. This explains a certain ambiguity in the representation (1.8). 
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7.2 Prom Jacobi to GMP 


In this section our goal is to estimate Pj{n) —Pj{oin) and qj{n) — qj{an) by means of the 
related distances dist((5“” J{E)) < oo and then apply Theorem 1.20. Therefore 

first of all we prove the following lemma. Note that the relation (7.5) below evidently 
implies a word-by-word counterpart of (1.3) in DKST, see Remark 1.7. 

Lemma 7.2. Let J be of the form (1.8). Then 

OO n 

y~]dist^((5~”'J5’^)+, J(an)+) < oo, an = '^ea{k) - pn. (7.5) 

n=0 fc=0 

Proof. We have 

n+k 

\b{k + n) - B{an - pk)\ < \eb{k + n)\ + Ci{E)\\ ^ eaij)\\, 

j=n+l 

where Ci{E) = sup^g^s/zs l|gradR(a)||. For p < 1, we have 

( n-\-k \ n-\-k n-\-k 

Ell E <EE E ik«(i)iP'^’(''' = EE E 

fc>l j=n+l J n>0k>lj=n+l k>ln>0j=n+l 


s E '‘P E *:|l««(i)ll" S E ■ E 

k>l j>l j>l k>l 

Making a similar estimation for \a{k + n)^ — A{an — pk)\ we obtain (7.5). □ 

Before to proceed we make the following important for us remark. 

Remark 7.3. Note that Az, see Lemma 5.6, is well define for all z G C \ cj(J). That 
is, in fact, we have a Hermitian analytic vector bundle in this domain. Its fundamental 
characteristic, the so-called curvature, is of the form A log(>r^, x^), see e.g. [8]. Being 
restricted on the real axis, it represents the Schwarzian derivative of r^{z). Our further 
considerations are based on estimations of related expressions and involve derivatives 
of exactly this level, see Lemma 7.4 below. We can conjecture that a certain Hermitian 
analytic vector bundle model, which generalized the model described in Section 2, is 
possible for operators of Killip-Simon class. Under more restrictive assumptions, when 
the absolutely continuous part of the spectral measure satisfies the Szego condition and 
positions of the eigenvalues outside E obey the Blaschke condition, such model does 
exist. This is the so-called scattering model for the given operator [28, 41, 27]. 

O 

Lemma 7.4. Let J G J{E). In the notations of Lemma 5.6, 

((J - J)xc, Xc) = sin((^(c) - V9(c)). (7.6) 
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Consequently, there exists C = C{a{J),c) < oo sueh that 

I sin((/9(c) — (/?(c))| < C'dist^(J+, J+) (7-7) 

and simultaneously for the derivatives 

|((^)(-)(c) - (^M(c)l < C'dist,(J+, J+) (7.8) 

for m = 1, 2,3 and rj > |6(c)|. 

Proof. We have 

((J - J)Xc, Xc) = {{J - c)Xc, Xc) - ((J - c)Xc, Xc). 


We simplify the first term 


1 


(e_ia( 0 ) sin (f + eo cos ip, Xc) = -((J — c)e_i, Xc) cos ip = sin cos ip. 

a(0) 

Thus, ((J — J)xc,Xc) = sincos (/9 — sin 99 cos and (7.6) is proved. 

The upper estimation in (5.15) in combination with (7.6), (5.14) implies 

|sin(^(c) -y(c))| < (7.9) 


Now, the vector —5—Xc in the functional model for J = J{a) corresponds to the normal 

Ikcll 

J_ ua 

fcpXc 

of this vector were estimated in (2.14). Thus 
1 


ized reproducing kernel h where Cc G D is such that 3(Cc) = c. The components 

II '' 


-||(J - J)xc|| < C'(S)dist^(J+, J+), |6(c)| < 7 ? < 1, 


Xp 


and (7.9) implies (7.7). 

To get (7.8) we differentiate (7.6) with respect to c 

COs(^(c) - </7(c))((^)'(c) - ^'(c)) = ((J - J)x', Xe) + ((J - J)Xe, (Xe)')- (7.10) 

Since sin(^(c) — ip{c)) was estimated from above, we have a uniform estimation for 
I cos((/9(c) — (/9(c)) I from below. Using (5.12), we evaluate x(.. Based on its explicit form 
and the estimation for ip'^, we obtain that ||x(.|| is also bounded by the distance from c 
to <j{J). Evidently, the coefficients of (xc)' also satisfies (2.14). Thus, 


|(^)'(c)-^'(c)| = 


||(J- J)x,||||x;,|| + ||(J- J)(xJ||||Xe 
|cOs((^(c) - (/9(c))| 


implies (7.8). Taking the second and third derivatives in (7.10), we obtain (7.8) for 
m = 2,3. □ 


42 



Corollary 7.5. If J is of the form (1.8) and A{n) =F ^{S ^JS^), then 


-Po(an)P < oo, an = '^ea{k) - fin. 


(7.11) 


n=0 A:=0 

Proof. By (7.7), (7.8) we can estimate the difference 

( 0 )^^^ „ ^ _ a(n)sinv9(ci) a(0) sin (^(ci) ° 


Po ’in)-poian) = 


J — J (Of),), 


7^'(ci) (^)'(ci) 

by means of dist((5“"'J5"')+, J{an)+). Due to (7.5), we have (7.11). □ 

Finishing the proof of Theorem 1.5. It remains to show that (1.8) imply (1.31)-(1.33). 
Similarly to (5.16), consider the ordered system of vectors 


e-i, x'ci, • • •, >!^cg,eo, ^ci) • • •) ^'cgFi- 
Let us point out that the orthogonalization of the system 

o o , o , ^ ° 

6—1, >ici i i ^Cg J CQ) (^Ci) ) • • • ) \^Cg ) ) Cl 

leads to the family see (2.6), where J = J{a). 


(7.12) 


(7.13) 


To evaluate the Gram-Schmidt matrix of the system (7.12) we use 

Cj - Cjn 

Therefore, 

, COs{(p{Cj) - (p{Cm)) ^ sin(v9(Cj) - piCm)) 

(^c,o =-:—- T [Cj) - 


[Cj - c 


m ) 


and 


V'-'j ) 

sin(v9(Cj) - p{Cm)) J 


+ 




ip {Cj^p (CjT^) 2 


sin((^(cj) - p[c^)) . 


yCj C'jYi) 


, j yim. 


Having uniform estimations from below for all Gram-Schmidt determinants of the 
system (7.13), from (7.7), (7.8), similarly to (7.11), we obtain 


^ \pp\n) - pj{an)\‘^ < oo, ^ \qp\n) - qj{an)\‘^ < oo, m = -1,0,1, j = 0,.. .,g. 

n=i) n=0 

This implies (1.31)-(1.33), in particular, ~< oo, j = 0, ...,g — l. 

□ 
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8 Appendix: one sided GMP matrices 


We describe interrelations between one sided GMP and Jacobi matrices given by, see 
(1.27), 

r_(2;) = ((J_ - 2:)“^e_i,e_i) = ((A_ - z)“^e_i,e_i) = f (gj) 

J X z 


Let us point out that both matrices have the same cyclic vector, and they are related 
by a common spectral measure da-. It is assumed that ci,..., do not belong to the 
(closed) support of this measure. 

Thus, as soon as A- is given we can construct J_ in the usual way, making basis of 
orthonormal polynomials. In the opposite direction we will construct an orthonormal 
basis of rational functions in . The matrix of the multiplication operator by the 
independent variable in this basis {Tk{x)}k<-i is A-. 


Definition 8.1. To the given da- we associate the orthonormal system 


T_i(x) r_2(x) 


T-g-l{x)_ = 


where L is the upper triangular matrix 


L = 


/-l) /- 2 ) 


^(- 9 - 1 ) 


^9 


L, 


> 0 , 


such that f Tk(x)Tj(x)da-(x) = 5kj■ 

In other words, if D is the Gram-Schmidt matrix of the given system 


( 8 . 2 ) 



1 

r-{cg) 

r_(ci) 


D = 

r-{cg) 

r'-{cg) 

r_(ci)-r_(cg) 

Cl—Cg 

r-{cj)-r-{ck) _ f da-{x) 




Cj-Ck J {X-Cj){x-Ck) 


_r_(ci) 

r- (ci)-r_ (cg) 

Cl—Cg 

^'.(ci) 



then L is defined via the upper-lower triangular factorization of D ^ 

L*DL = / or D-^ = LL*. (8.3) 

Lemma 8.2. The multiplication operator with respect to the orthonormal system (8.2) 
represents the B-bloek of a GMP matrix (see (1.22), (1.23 )), that is, for the matrix B-i 
given by 

•= j XTj{x)Tk{x)da-{x) (8.4) 
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its lower triangular part (including the main diagonal) is of the form 


B+ = (m(£)*)+ + C, (8.5) 

where 

f T-i(x)xda-(x) 0 

f T-2(x)xda-(x) ^ Cg 

£ := ^0 ,m= , C = 

£(-9-1) j T-g-i{x)xda-{x) Cl 

Proof. Note that 

1 ^ ••• ^ = [^ -1 ••• -1] + 1 ^ ••• ^ C. (8.6) 

We substitute this and ( 8 . 2 ) in (8.4). Since constants are orthogonal to Tk{x) for all k, 
except for A: = —1, we obtain 

B = m £(-2) ... + ho [0 * ... *] + L*DCL. 

By (8.3) L*DCL = p-^CL. Since {L-^CL)+ = C, we have (8.5). □ 

Now, we consider the system (8.2) as a cyclie subspace for the multiplication by 

9 ^ 

A(x) = Aqx + Co + ^ Xj >0, j = 0,...,g, (8.7) 

k=l 

in . We inductively define the orthonormal system 

T_m{g+i)-k-i{x) = + ..., m>l, 0 <k < g. ( 8 . 8 ) 

Lemma 8.3. The multiplication by A(x) with respect to the system ( 8 . 8 ) is a (2g + 3)- 
diagonal matrix, or a {g + 1)-block diagonal Jacobi matrix 

'tu-i Dhi 

g_ = t)_i tu_2 0^2 , ( 8 . 9 ) 

with upper-triangular having positive diagonal entries. 

£(’T‘) 

Proof. The relation Q-Cm = ^m-g-i + ■ • • follows from the definition (8.8). That 

^-m+g 

is, dm is upper-triangular. The fact that the operator is self-adjoint implies that all dm 
are lower-triangular matrices. 

Note that corresponding to this Q- matrix measure d'P,-{y) is of the form (5.5). □ 
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As it was claimed, we define A- = A-(J-) as the matrix of the operator multipli¬ 
cation by X in the basis (8.8). Evidently, t/_ and A_ commute. 

Theorem 8.4. In the given eonstruetion A_ = A-{J^) is a one sided GMP matrix. 


Proof. In Lemma 8.2 it was checked that the block i?_i has the required structure. We 
claim that the same is related to the block A_i. Indeed, the function A(x)r_i(x) = 
A(x)£q is a linear combination of x and functions from the chosen cyclic subspace, 
see (8.7). That is, T-g- 2 {x) is of the same form. Thus, for all j = —1,..., —g — 1, by 
(8.6 ), XTj{x) is a linear combination of {Tj(x)}~fZ^. In other words A_i is of the form 


A-i = So 


Po 


pi 


= Soip-iY 


( 8 . 10 ) 


where 


{p-iY := j T_g_ 2 {x)x[T-i{z) T- 2 {x) ... T-g-i{x)] da{x) 


= / r_ 3 _ 2 (x)([x -1 ... -l] + 


1 


pi~‘A 

'0 ^0 ■ ■ ■ ^0 


T_g_2{x)xda{x) = 


Cl —X 

(-1) d-2) 


C)Lda{x) 


0 


^0 


^0 


\ d-l) 4-3-2) 

MIq ^g+1 


It remains to use the commutant relation Q-A- = A-Q-. We have 


Due to (8.10), we get 

A -2 = (o_ 2 )oo<^o(P-i)*ol}, that is, A _2 = Soip- 2 )*, {P- 2 )* = (0-2)oo(P-i)*ol}. 
Generally, 

Ak = So{pk)*, where (pk)* = (Ofc)oo(A+i)*0^+i. (8.11) 

For B-blocks we have 

Bk — (ffc+iT^fc+i “1“ . 

If we assume 

Bk+l = Qk+l{Pk+lT + 

where Mk+i is upper-triangular, which main diagonal is C, then 

Bk = t)k+iqk+i{pk+i)*'o^li+tVkSo{pk+i)*'o^li+Mk = {t)k+iqk+i+tOkSo){pk+i y^kli+^k 
where Mk is also upper-triangular. From this relation and (8.11) we get 

Bk = qkiPkT +Mk, 

where up to the first component the vector ijk has the form \ ^ {t>k+iqk+i TtOfedo) and 
by definition Mk preserves the structure of its main diagonal for all k. Note that this is 
possible, since the term = (lo(pfc-i-i)*tDfc-i-iO^_^;^ has a non vanishing entry 

on the main diagonal only in the first component. 

□ 
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